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Microwave observations of spinning dust emission in NGC 6946^ 
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ABSTRACT 

We report new cm- wave measurements at five frequencies between 15 and 18 GHz of the con- 
tinuum emission from the reportedly anomalous "region 4" of the nearby galaxy NGC 6946. 
We find that the emission in this frequency range is significantly in excess of that measured 
at 8.5 GHz, but has a spectrum from 15-18 GHz consistent with optically thin free-free emis- 
sion from a compact Hll region. In combination with previously published data we fit four 
emission models containing different continuum components using the Bayesian spectrum 
analysis package radio spec. These fits show that, in combination with data at other fre- 
quencies, a model with a spinning dust component is slightly preferred to those that possess 
better-established emission mechanisms. 

Key words: radiation mechanisms: general — galaxies: individual(NGC 6946) 



1 INTRODUCTION 

The complete characterization of microwave emission from spin- 
ning dust grains is a key question in both astrophysics and cosmol- 
ogy as it probes a region of the electromagnetic spectrum where a 
number of different astrophysical disciplines overlap: it is impor- 
tant for CMB observations in order to correctly characterise the 
contaminating foreground emission (Gold et al. 2010); for star and 
planetary formation it is important because it potentially probes 
a regime of grain sizes that is not otherwise easily observable 
(Rafikov 2000). 

Although a number of objects have now been found to ex- 
hibit anomalous microwave emission, attributed to spinning dust, 
it is still unclear what differentiates those objects from the many 
other seemingly similar targets that do not show the excess. In or- 
der to investigate this question a number of Galactic observations 
have been made towards known star formation regions (see e.g. 
AMI Consortium: Scaife et al. 2010 and references therein; Casas- 
sus et al. 2008; Tibbs et al. 2009). In addition Murphy et al. (2010; 
hereinafter MIO) made the first extra-galactic search for anomalous 
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microwave emission within the star formation regions of the nearby 
galaxy NGC 6946 using the Caltech Continuum Back-end on the 
Green Bank Telescope. MIO found a significantly anomalous spec- 
trum in only one of 10 star-forming regions: extra-nuclear region 4 
(hereinafter NGC 6946-E4). The excess of emission was seen be- 
tween 27-38 GHz relative to the continuum emission at 8.5 GHz 
measured using combined Effelsberg 100 m Telescope and VLA 
observations (Beck 2007). 

In this letter we present follow-up observations of NGC 6946- 
E4 at frequencies from 15-18 GHz using the Arcminute Mi- 
crokelvin Imager (AMI) Large Array (LA). In Section 2 we present 
the details of these observations, in Section 3 we present the results 
of the AMI-LA observations and a comparison with other radio 
data, and in Section 4 we discuss the implications of these results 
and form our conclusions. In what follows we use the convention: 
»S oc v", where S is flux density, v is frequency and a is the spectral 
index. All errors are quoted to 1 CJ. 



2 OBSERVATIONS 

The AMI-LA consists of eight 13 m antennas and is sited at Lord's 
Bridge, Cambridge (AMI Consortium: Zwart et al. 2008). The 
telescope observes in the band 13.9-18.2 GHz with cryogenically 
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cooled NRAO indium-phosphide front-end amphfiers. The overall 
system temperature is approximately 25 K. Amplification, equal- 
ization, path compensation and automatic gain control are applied 
to the IF signal. The back-end has an analogue lag correlator with 
16 independent correlations formed at double the Nyquist rate for 
each baseline using path delays spaced by 26 mm. From these, eight 
complex visibilities are formed in 0.75 GHz bandwidth channels. 
In what follows, the three lowest frequency channels (1, 2 & 3; 
13.9-14.6 GHz) are not used due to interference from geostation- 
ary satellites. 

NGC 6946-E4 (J 20^34"^ 19.^7 + 60°10'08:'7) was observed 
by the AMI-LA in a single 12 hour synthesis. Data reduction was 
performed using the local software tool REDUCE. This applies both 
automatic and manual flags for interference, shadowing and hard- 
ware errors, phase and amplitude calibrations and then Fourier 
transforms the correlator data to synthesize frequency channels be- 
fore output to disk in uv FITS format suitable for imaging in MPS. 
Flux calibration was performed using short observations of 3C48 
at the beginning of the run and 3C286 at the end of the run. We 
assumed I-i-Q flux densities for these sources in the AMI-LA chan- 
nels consistent with the frequency dependent model of Baars et al. 
(1977), ~ 1.64 and 3.48 Jy, respectively, at 15 GHz. As Baars et al. 
measure I and AMI-LA measures I-i-Q, these flux densities include 
corrections for the polarisation of the calibrator sources derived by 
interpolating from VLA 5, 8 and 22 GHz observations. A correc- 
tion is also made for the changing intervening airmass over the ob- 
servation. From cross-calibration of 3C48 and 3C286, we find the 
flux calibration is accurate to better than 5 per cent. The phase was 
calibrated using interleaved observations of J 203 1+5455, selected 
from the Jodrell Bank VLA Survey (JVAS; Patnaik et al. 1992). 
After calibration, the phase is generally stable to 5° for channels 
4-7, and 10° for channel 8. The FWHM of the primary beam of the 
AMI-LA is?^6'atl6GHz. 

Reduced data were imaged using the AlPS data package. 
CLEAN deconvolution was performed using the task IMAGR. 
CLEAN deconvolution maps were made from both the combined 
channel set and for individual channels. 



3 RESULTS 

3.1 AMI Large Array Data 

The combined channel map from the AMI-LA observation of 
NGC 6946-E4 is shown in Fig. 1. E4 is located at the phase 
centre with extra-nuclear region 8 (J20^34"^32.^52 + 60°10'22f'0; 
NGC 6946-E8) to the east and the galactic nucleus of NGC 6946 
(J20h34"^52.^34 + 60°09'14f'2) further to the east, outside the 
primary beam FWHM. Flux densities were extracted for both 
NGC 6946-E4 and NGC 6946-E8. Since NGC 6946-E8 also lies 
within the FWHM of the primary beam we used the recovered 
flux densities for this source as a check on the absolute calibra- 
tion for this field. The flux densities measured for NGC 6946-E4 
and NGC 6946-E8 are listed in Table 1. 



3.2 NCC6946-E4: Comparison with other radio data 

We compared the AMI-LA data with that at 8.5 GHz from the Ef- 
felsberg 100 m and VLA telescopes (Beck 2007). These data in 
their original form, see Fig. 1(b), constitutes a total power measure- 
ment of the region at a resolution of 15''. We Fourier transformed 
these combined data and sampled them at uv positions identical to 



those of the AMI-LA 15.0 GHz data (channel 4) in order to com- 
pare consistent angular scales. The resulting uv data set was then 
mapped and cleaned in the same way as the AMI-LA data. Since 
channel 4 of the AMI-LA recovers a large amount of the extended 
structure around NGC 6946-E4 the flux density for this source from 
the sampled 8.5 GHz map is very similar to the unsampled value, 
see Table 1. Flux densities were extracted from both the AMI-LA 
and sampled Effelsberg 8.5 GHz maps using the FITFLUX software 
(Green 2007; AMI Consortium: Scaife et al. 2009). This method 
calculates flux densities by removing a tilted plane fitted to the lo- 
cal background and integrating the remaining flux. We do this by 
drawing a polygon around the source and fitting a tilted plane to the 
pixels around the edge of the polygon. Where an edge of the poly- 
gon crosses a region confused by another source the background 
is subjective and we omit this edge from the fitting. An example of 
where this might be appropriate is shown in Fig. 1 (b). Since the ex- 
tracted flux density is dependent to some degree on the background 
emission we repeat this process using five irregular polygons, vary- 
ing each slightly in shape. The final flux density is the average of 
that extracted from these five polygons. 

Errors on the flux densities were calculated as CTy = 

^J (0.055'v)^ + 0'^s,v + <^fit' where <7rms,v is the rms noise outside 
the primary beam on each channel map and Gf^i is the standard de- 
viation of the fluxes measured in the five polygonal apertures. The 
errors are dominated by CJfit, which is large due to the complicated 
background emission in this crowded field. 

Since the uv coverage of the AMI-LA varies across the fre- 
quency channels we quantified the amount of flux lost in each chan- 
nel relative to channel 4 by sampling the total power 8.5 GHz map 
to match the uv coverage in each of the AMI-LA frequency chan- 
nels. This showed that flux loss was negligible in channels 4-6 and 
notable (^ 5%) in channels 7 & 8 only. We corrected for this loss in 
the measured AMI-LA flux densities with corrected values shown 
in parentheses in Table 1 . In the same manner we checked for cor- 
rections to the flux densities of NGC 6946-E8. In this case the cor- 
rections were negligible (< 2%) and no corrections were made. We 
found the 15-18 GHz data for NGC 6946-E8 to be consistent with 
the spectrum for this object presented in MIO. The flux density of 
NGC 6946-E4 across the AMI-LA band is in excess of the 8.5 GHz 
flux density by approximately 10 cr. 

The collected observations of both region E4 and E8 were 
analysed using an updated version of the radio spec package^ 
(Nikolic 2009). This software tool calculates the Bayesian poste- 
rior distribution of the parameters of a model for the radio spec- 
trum and, furthermore, the Bayesian evidence. The implementa- 
tion of these calculations is based on the nested sampling algo- 
rithm by Skilling (2004). For the analysis of the data in this pa- 
per, we used several different models which each consist of a num- 
ber of components, each with physically parameterized properties. 
Two components that are present in all models are a synchrotron 
component, which we parametrise in terms of the supernova rate 
within the beam, and an un-absorbed free-free component, which 
we parametrise in terms of the star-formation rate within the beam. 
The conversions from these physical parameters to radio luminosi- 
ties are made according to the formulae given by Condon (1992). 

In order to explain the excess of emission at cm wavelengths, it 
is necessary to introduce another component that contributes to the 

^ The complete code and data used for this spectrum anal- 
ysis are available for public download under GPL license: 
http : / / www . mrao . cam . ac . uk/ ~bii204/ galevol/ speca/ sdgals . html. 
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Figure 1. NGC 6946-E4. (a) Greyscale and contours at —3,3,6, 12 cr etc are from the AMI-LA where arms = 25 /iJy beam~^, greyscale units are mJy beam^^ . 
Data are not corrected for the primary beam response. The AMI-LA synthesized beam, 29.2" x 25.2", is shown as a filled ellipse in the bottom left corner, (b) 
Contours are from the AMI-LA and greyscale at 8.5 GHz is from Effelsberg-VLA measurements, greyscale units are /iJy beam^ (Beck 2007). The positions 
of NGC 6946-E4 (at the phase centre) and NGC 6946-E8 (to the east) are marked with crosses, and a line is drawn to indicate a non-fitted polygon edge (see 
text for details). The AMI-LA primary beam is shown as a black circle in both images. 



Table 1. Flux densities. 



AMI-LA Channel Number 



Source 8.5 GHz 4 5 6 7 8 

(mJy) (mJy) (mJy) (mJy) (mJy) (mJy) 

[2] [3] [4] [5] [6] [7] 

NGC6946-E4 1.94 ±0.29 3. 17 ±0.36 3. 16 ±0.33 3.00 ±0.28 3. 15 ±0.43 2.78 ±0.34 

(3.31 ±0.43) (3.09±0.35) 

NGC6946-E8 2.11 ±0.21 1.90 ±0.27 1.51 ±0.27 1.40 ±0.23 1.51 ±0.34 1.42 ±0.19 

Column [1] contains the source name, [2] the wv-sampled 8.5 GHz flux density for that source, [3]-[7] contain the AMI-LA flux densities for 
each source from AMI-LA channels 4-8, with flux loss corrected values in parentheses.. 



emission. As outlined above we have considered two options: emis- 
sion due to spinning dust, and absorbed free-free emission from 
a compact H II region. For the spinning dust emission component 
we used the warm ionized medium model described by Draine & 
Lazarian (1998)^. The only degree of freedom in this model is the 
overall amplitude, which we parametrised in terms of the total mass 
of gas carrying the spinning dust. 

Our model for the absorbed free-free emission from H II is 
again parametrised in terms of the star-formation rate in the region, 
but it now also has the filling factor, i.e. the fraction of the area 
covered by the telescope beam that the H ii region subtends, as a 
free parameter. This area is used to compute the electron opacity 
as a function of frequency and to correct the unabsorbed free-free 
model for the effects of absorption. 



4 DISCUSSION AND CONCLUSIONS 

Considered on their own, the AMI-LA data (after correction for 
flux loss) have a spectral index of Of ami = —0.11 ±0.77. Although 
this value is consistent with optically thin free-free emission, the 

^ ftp : / /ftp . astro . princeton . edu/ draine/ dust/ 
spin/ emit4 . j nu . wim_a 



error on this estimate is large and we cannot rule out other mecha- 
nisms. Since the spectral index between the Effelsberg-VLA mea- 
surement at 8.5 GHz and the AMI band is rising (a^^ = 0.67 ± 
0.08, see Fig. 2), we need to consider the possibility that region E4 
contains one or more compact Hll regions with their opacity reach- 
ing unity at approximately 12 GHz. Such an opacity would require 
an emission measure of ~ 5 x 10^ pc cm~^, assuming Tq = 10^ K, 
and would be appropriate for a compact Hll region. 

We therefore examine two alternative hypotheses for the emis- 
sion from region E4. The first is that the emission is due to the usual 
diffuse synchrotron and free-free mechanisms associated with star- 
formation, with an additional high-opacity free-free component 
(Hypothesis 1; HI). The second hypothesis is that there is a spin- 
ning dust component rather than high-opacity free-free (H2). A 
summary of how well these two hypotheses fit the observed data 
is shown in the form of fan-diagrams in Fig. 2. As can be seen from 
this figure, neither of the hypotheses can be ruled out, although the 
spinning-dust appears to somewhat better match the data. This is 
also confirmed by a simple comparison of the models: assuming 
flat priors and no a priori difference between the models, the log- 
arithmic Bayes factor is 3.7 ±0.3 in favour of the spinning dust 
model. From the Jeffreys' scale of evidence (Jeffreys 1961; Kass 
& Raftery 1995; see e.g. Efron & Gous 2001 for further discussion 
of this scale) this would indicate a weak positive preference for 
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Figure 2. The observed radio spectrum of region E4 of NGC 6946 (points and error bars) with the fan-diagram of two model fits to these data: on the left is the 
model with a highly absorbed free-free emission region (HI) and on the right is the model with spinning dust emission (H2). Low frequency data are taken 
from MIO, scaled to the wv-sampled flux density at 8.5 GHz, with the exception of points between 15 and 18 GHz which are from the AMI-LA (this work). 
The assumed error is 10% unless stated otherwise in Table 1. The colour scale indicates the evidence contribution as a function of frequency and flux density, 
for details see Nikolic 2009. 
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Figure 3. The observed radio spectrum and fan-diagrams for region E8 of NGC 6946, with data and errors as in Fig. 2. The model on the left only consists of a 
synchrotron and un-absorbed free-free components (HI) while the model on the right also has a spinning dust component (H2). Low frequency data are taken 
from MIO, scaled to the wv-sampled flux density at 8.5 GHz, with the exception of points between 15 and 18 GHz which are from the AMI-LA (this work). 
The assumed error is 10% unless stated otherwise in Table 1. Colour scale as above. 
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Figure 4. Marginalised distribution of inferred gas mass that bears the spinning dust: the left panel shows extra-nuclear region 4 and the right panel extra- 
nuclear region 8. 
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Table 2. Derived model parameters and errors 



Model 


SNe rate 


O^sync 






^gas 


/ 




logio(yr"^) 




logio(yr~^) 


logio(yr~^) 


logio(M0) 


logio(str) 




[2] 


[3] 


[4] 


[5] 


[6] 


[7] 


Prior 


(-5,-1) 


(-1.0,-0.5) 


(-3,-1) 


(-3,0) 


(5,9) 


(-6,-3) 


NGC 6946-E4 (HI) 


-3.84±0.11 


-0.74 ±0.14 


-0.93 ±0.07 


-1.80±0.48 




-4.97±0.18 


NGC 6946-E4 (H2) 


-3.98±0.16 


-0.71 ±0.15 




-1.40 ±0.36 


8.12±0.09 




NGC 6946-E8 (HI) 


-3.65 ±0.04 


-0.66±0.13 




-1.70±0.41 






NGC 6946-E8 (H2) 


-3.64 ±0.04 


-0.65±0.12 




-1.75±0.41 


6.19±0.66 





the spinning dust model above the free-free model. The maximum 
likelihood parameters for these models are listed in Table 2. 

For comparison we have also carried out a similar analysis 
on the region E8, shown in Fig 3. In this case the two hypotheses 
are a simple diffuse synchrotron plus free-free model (HI) versus 
the same model with an additional spinning dust component (H2). 
In this case the logarithm of the Bayesian evidence ratio is 0.5 ± 
0.3 in favour of the simpler model without the spinning dust. A 
ratio of this size indicates no perceptible difference between the 
two models. 

In region E4 where a spinning dust model is the preferred hy- 
pothesis, we can marginalise the posterior distribution of the model 
parameters to obtain an estimate of the gas mass containing the 
spinning dust, shown as a histogram in the left panel of Fig. 4 and 
in numerical form in Table 2, i.e., iQ^ l^^ l M©. In region E8 spin- 
ning dust is not the preferred hypothesis but proceeding with this 
hypothesis anyway, we find an upper limit for the mass of the gas 
bearing the spinning dust, which is around 10^ -^Mq. We draw a 
conclusion that if the conditions in E4 and E8 are similar, then the 
mass of any gas bearing spinning dust in E8 must be at least a factor 
of five smaller than in E4. 

Region E4 is located on the dense rim of a "remarkable" Hi 
hole (Boomsma et al. 2008) within NGC 6946. Such an association 
may be relevent to the differentiation of this star formation region 
from the eight others found to exhibit no anomalous emission by 
MIO. The hole is remarkable for a number of reasons, notably the 
almost unbroken symmetry of its dense Hi rim, unusual in so large 
a hole, and the small scale high velocity gas complexes observed in 
connection with it. 

As mentioned above, the spinning dust model is preferred by 
the evidence calculation for NGC 6946-E4, but not at a very high 
level. Definitive confirmation of the nature of the emission requires 
measurements at frequencies above 50 GHz where the spinning 
dust and compact Hll region models have significantly different 
behaviour. For example. Fig. 2 shows that at 100 GHz the differ- 
ence between these two models should be at least a factor of two in 
brightness. 

In the sub-mm there are data available from SCUBA at 
850 jUm for this region by Di Francesco et al. 2008, which might 
be used to constrain the mass of NGC 6946-E4 and place con- 
straints on the frequency at which the optical depth reaches unity. 
If a compact Hll region is present, with T = 1 at V > 8.5 GHz 
it should have a correspondingly high dust mass. From analysis 
of the SCUBA data we obtained a flux estimate for region E4 of 
5850 = 11 ± 14mJybeam~^. However, the errors on this estimate 
are too high to allow any useful constraint on the properties of the 
thermal dust emission or to calculate a reliable dust mass estimate. 

From the existing data the possibility of a spinning dust com- 
ponent in this region cannot be ruled out, but the evidence is not yet 



definitive. Further observations of this object at frequencies cover- 
ing the higher frequency minimum between spinning dust emission 
and thermal dust emission 90 GHz) would be most useful. 
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